Kondo effect in a quantum dot coupled to ferromagnetic leads: A numerical 

renormalization group analysis 
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We investigate the effects of spin-polarized leads on the Kondo physics of a quantum dot using 
the numerical renormalization group method. Our study demonstrates in an unambiguous way that 
the Kondo effect is not necessarily suppressed by the lead polarization: While the Kondo effect is 
quenched for the asymmetric Anderson model, it survives even for finite polarizations in the regime 
where charge fluctuations are negligible. We propose the linear tunneling magnetoresistance as an 
experimental signature of these behaviors. We also report on the influence of spin-flip processes. 
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Introduction. — Magnetic impurities embedded in 
metallic hosts cause anomalous resonant scattering of 
conduction band electrons. At the same time, the local- 
ized magnetic moments are screened at low temperature 
by the itinerant electron spins. This is the celebrated 
Kondo effect 0, which has been recently revived in 
mesoscopic physics Q • Ever since the theoretical predic- 
tions Pi LJ and the experimental demonstrations |a] , the 
Kondo effect in phase-coherent systems such as quantum 
dots (QD's) has stimulated great interest in this field. 
The remarkable success behind this is the fine tunability 
of the parameter space (impurity level and hybridization 
couplings). The controlled manipulation in mesoscopic 
systems has not only allowed to test various aspects of 
the Kondo effect, which is a hard task in bulk solids, but 
also has posed further exciting questions. For example, 
when the spin- degeneracy of the impurity level is lifted 
by an external magnetic field, the Kondo peak in the 
density of states (DOS) of the dot is expected to split |fj. 
However, new experiments Q and theoretical studies 
suggest that the situation is more subtle. 

A flood of very recent works H ES [H E H H 
has introduced another interesting issue, namely, 
how the Kondo physics is affected when the continuum 
electrons themselves are allowed to form spin- dependent 
bands. The motivation for this research stems from the 
successful field of spintronics 01- I n particular, a change 
has been detected in the resistivity of a Kondo alloy due 
to spin-polarized currents |isj . Furthermore, it is already 
possible to attach ferromagnetic leads to a carbon nan- 
otube 0| , and a carbon-nanotube QD has been shown to 
display Kondo physics below an unusually high temper- 
ature .20]. fn addition, a QD coupled to ferromagnetic 
electrodes has been proposed as a promising candidate 
for spin injection devices, but so far, analyzed only in 
the Coulomb blockade regime . The present work pro- 
vides precise theoretical predictions in a wider region of 
the parameter space including the strong coupling limit. 

The aim of our work is twofold. First, based upon a 
numerical renormalization group (NRG) calculation, we 



investigate the influence of ferromagnetic electrodes and 
the relative orientations of their magnetizations on the 
equilibrium properties of a QD with and without intrin- 
sic spin flip processes. To the best of our knowledge, 
this is the first study of the model that sweeps across the 
different regimes (i.e., Kondo, mixed- valence, and empty- 
orbital), thoroughly analyzed with the assessment of lo- 
cal DOS, linear conductance, and tunneling magnetore- 
sistance (TMR). Second, we resolve a controversy lately 
raised in the literature with regard to whether a spin- 
dependent renormalization of the impurity level induced 
by the spin-polarized leads will split the Kondo peak 
when the magnetic moments of both leads are aligned. In 
Ref. ||, an equation-of- motion (EOM) method plus an 
ansatz for the interacting self-energy |22| were employed 
and it was suggested that the splitting, 5, is absent. In a 
later work , the scaling arguments (together with the 
EOM method) were used to find that S is nonzero. In 
Ref. [Ioj |. using a similar approach, a splitting was pre- 
dicted only in the mean-field peaks. In a more recent 
work |l6j , they made use of a noncrossing approximation 
(NCA) and obtained 6 ^ 0. In Refs. [HO, on the other 
hand, the slave-boson mean-field theory (SBMFT) was 
utilized to study the zero-temperature properties and no 
splitting was observed. The answer to the controversy is, 
thus, elusive because each approximation method men- 
tioned above has certain drawbacks of its own. Below, 
according to a NRG calculation, which has been known to 
provide very accurate results for impurity problems |'2.'"{l| , 
we resolve clearly the controversy. 

Our main result is summarized in Fig. EJa). We find 
that in the presence of electron- hole symmetry the Kondo 
peak at the Fermi level remains unsplit even at finite po- 
larizations and the linear conductance achieves the uni- 
tary limit. This remains true as long as only spin fluctua- 
tions are present in the QD. On the contrary, when charge 
fluctuations start to play a role (as in the asymmetric case 
of the Anderson model), the Kondo peak shows a visible 
splitting and the conductance is then suppressed. 

Model. — We consider an ultrasmall tunnel junction 
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FIG. 1: Left: Local DOS of the QD for the symmetric 
Anderson model, (a) Af(o>), (b) and (c) A(oj) for 

e d = -U/2 = -0.1D (D is the bandwidth). Right: DOS for 
the asymmetric case, (d) A-[{uj), (e) Ai{uj), and (f) A(uj) for 
e d = -0.1D and f7 = 0.45.D. In all cases T = 0.02D. 

comprising a QD coupled to two ferromagnetic leads. 
We assume that the QD consists of an energy level with 
an unpaired spin- 1/2 electron and a charging energy 
U. This way, the QD is equivalent to an Anderson- 
type impurity with single-particle energy Ed,a for spin 
CT = {1\ 1} |ia L| Notice that £d,<j includes the Zeeman 
energy Az = £d,T — e d,i of an external magnetic field. In 
what follows, we set Az = in order to unmask possible 
spin-dependent renormalizations of the bare energy level 
purely due to coupling with the leads |24| . On the other- 
hand, we shall include in a phenomenological way inter- 
nal spin- flip scattering processes with rates 2R/h |25|. 
Tunneling of electrons from the QD to the leads (reser- 
voirs) a — {L, R} is described by the hopping integral 
V a ka- The resulting Hamiltonian is given by: 

J4? = ^ Ed,a^<T + Uh^hi + R (d^d^ + H.c.J 

a 

+ £<xk*cl k(7 C a ko + ^{Votkacl^da + H.C.) , (1) 
aha aka 

where c^ fco . (c Q fc CT ) is the creation (annihilation) operator 



for an electron with wave vector k and spin in the elec- 
trode a. The QD occupation number is h a = d\d a [d^ 
(d CT ) creates (annihilates) an electron in the dot]. 

For definiteness, we shall take identical leads with 
chemical potentials [ip = = Ep and symmetric cou- 
plings (i.e., equal tunnel barriers). Ferromagnetism on 
the leads may be represented either by a spin-dependent 
DOS /9 aC r(w) or by spin-dependent tunneling amplitudes 
Vaka ■ We choose the latter for convenience but both pic- 
tures are formally equivalent provided we are interested 
just in the transport properties of the junction. In any 
case, the overall effect results in a spin-dependent hy- 
bridization parameter T aa (uj) = ir J^k \Vaka\ 2 3(u — £ak)- 
(We neglect proximity effects such as stray fields induced 
in the QD). A further simplification consists of neglect- 
ing the energy dependence of T aa (uj), evaluating it at 
id = Ep (wide band limit). In the following, we choose 
Ep = as the origin of energies. Therefore, one could 
define the spin polarization (close to the Fermi energy) 
at each lead as p a = (r Q | — T a i)/(T a ^ + T a i) with 
— 1 < p a < 1. We consider parallel (P) and antipar- 
allel (AP) magnetizations of the two leads. In the P 
case (p L = p R = p), we have T LT = r_R T = (1 +p)r /2 
and T Ll = T RL = (1 - p)T /2, where T = T qT + T a ±, 
whereas the AP case (pl = ~Pr = p) yields T p^ = T R i = 

(i + P )r /2andr^ = r RT = (i-p)r /2. 

In order to apply the NRG technique more efficiently, 
we map Eq. JQl onto an effective model with a single 
lead. This is achieved by means of a canonical transfor- 
mation 0, |2|| . In the P configuration it reads (we omit 
the index k for simplifying notation): 

CLa = + b a )/V2 , c Ra = (a a - b a )/V2 . (2) 

For AP polarizations one uses the relations: 

cl T = (V a ai + Vib^/V , c m = (V,a T - K a 6;)/V , (3a) 
CL l = (V l a l + V a b l )/V, c Rl = (V aai - VibJ/V, (3b) 

where V a = Vp-\ = Vri (K: = Vli = V R -\) is the tun- 
neling amplitude for majority (minority) spins and V = 
V\ y a\ 2 + \Vi\ 2 - Substituting Eqs. ©-© into Eq. Q one 
can show that (i) the QD electron decouples from the fee- 
operators and hybridizes only with the quasiparticles de- 
scribed by the a a operators; and (ii) the effective dot-lead 
couplings are renormalized for both configurations. 

We use the NRG method to obtain the local DOS 
A a (u>) = A>{lj) - A<(lu) on the QD (i.e., at the im- 
purity site). The particle and hole branches, A> (ui) and 
A<(lj), are defined by 

A>/< H = ±£| (m| dl/d a |0) | 2 5{u> T E m ± E ) , (4) 

1 a 

where |0) is the many-body ground state (electrodes + 
QD) and |m) is an excited state (E and E m are the cor- 
responding energies). Notice that all the physics (correla- 
tions, dependence on the gate voltage, etc.) is contained 
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FIG. 2: (a) Splitting S of the Kondo peak as a function of Ed 
for p = 0.25 and U = OAD. (b) S versus p for e d = — O.LD and 
U = 0A5D. (c) TMR versus p for e d = -0.1D and (7 = OAD. 
(d) TMR as a function of e d for p = 0.25 and (7 = OAD. In 
all cases T = 0.02D. 



in A a {uS). The linear conductance (normalized to e 2 /h) 
of the junction at zero temperature is obtained from the 
impurity spectral density function at the Fermi level |27| , 
.9 = E* 2T L*T Ra A a {0)/{r La + T Ra ). As it is the case 
in usual Kondo experiments on QD's, we shall present 
results for large interaction (U 3> F). 

Results. — We first address the issue whether the Kondo 
peak in local DOS splits or not. (Until later, we put 
R = 0.) Figure [2 shows A^(ui), A^uj), and A(uS) = 

(u>)+ A I (uj) for different values of the lead polarization 
p in the P configuration. [The AP case is less interesting 
as both spin orientations are equally coupled after the 
transformation given by Eq. |jSJ] . Left panels correspond 
to the symmetric Anderson model (i.e., Ed = —U/2). For 
p = 0, in addition to two (symmetric) mean-field peaks 
at ui — Ed and ui = Ed + U, A(u) shows a peak at ui — 
[see Fig. Hfc)], which is responsible for the observed 
zero-bias anomaly. As p increases, the spectral peak of 
A^(uj) [Ai(u>)] at the Fermi energy increases (decreases) 
[see insets of Fig. njb) and (c)]. Remarkably, however, 
the central peaks of both A^(uj) and A^(tu) are pinned at 
the Fermi level; in particular, the Kondo peak in A(uj) 
does not split. Experimentally, one would see a perfect 
transparency of the junction (see below). Right panels 
of Fig. H show the same functions for the asymmetric 
case (Ed 7^ — U/2), where charge fluctuations are allowed 
to certain extent. As p increases, A^ and A± shift in 
opposite directions [see Fig.^d) and (e)] and the Kondo 
peak in A(lo) splits into two [Fig.^f)]. As a result, the 
Kondo effect is suppressed. We have checked as well that 
both mean-field peaks are shifted in opposite directions, 
though their strengths differ and the splitting cannot be 
resolved in Fig. ^f). 
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FIG. 3: The effects of spin-flip processes R for the symmetric 
(left panels with Sd = —U/2 = — 0.1D) and asymmetric (right 
panels with Ed = —0.1D and U = 0.4D) Anderson model. 
Local DOS for (a,d) p = and (b,e) p = 0.25. (c,f) TMR 
versus R for p = 0.25. In all cases Fq = 0.02D. 



The finite splitting for the asymmetric Anderson model 
may be understood in terms of simple scaling argu- 
ments [l3j |: Because the hybridization for up spins is 
larger than for down spins (ro,| > r Q ^), the renormal- 
ization of the bare level Ed is spin-dependent; the t(l) - 
electron lowers (raises) its energy. Then, the coupling 
acts as an effective magnetic field, leading to a finite 
5 28]. Yet, the perturbative nature of a poor man's 
scaling cannot describe the fixed point in the strong 
coupling regime. In particular, such simple scaling ar- 
guments cannot account properly for the particle-hole 
symmetry in the symmetric Anderson model and always 
predict S ^ 0. For the symmetric Anderson model, it 
is important to notice that the particle-hole symmetry 
quenches charge fluctuations completely for both spins 
(( n t) = ( n l) = 1/2) at any \p\ < 1, and the real part 
of the self-energy (at Ep ) is zero. This means that al- 
though the binding energy of the singlet state (the Kondo 
temperature Tk) diminishes withp, the quasiparticle life- 
time is still infinite and the Fermi liquid picture is valid. 
Therefore, the results in Fig. ^c) are consistent with 
SBMFT, which describes the Kondo peak when spin fluc- 
tuations prevail. Likewise, the results in Fig. ^f) are in 
agreement with EOM and NCA models, which support 
charge fluctuations to some degree. Of course, the NRG 
method can encompass the whole parameter range. 

To illustrate our conclusions, we measure the splitting 
S of the Kondo peak as a function of Ed (experimen- 
tally this is controlled by a gate voltage) with U fixed 
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[see Fig. Ufa)]. The splitting S increases roughly lin- 
early from zero as moving away from the symmetric point 
Ed = —U/2. Notice also that well away from e d = — U /2, 
(5 is linear in the lead polarization [see Fig. 12(b)], con- 
firming the prediction relying upon scaling arguments. 

We now turn to the tunneling magnetoresistance of the 
system, defined by TMR = (g p — g 
show that g p = nT [(1 +p)A|(0) - 



AP )/g AP - One can 
{l-p)A l (0)] and 



g AP = (1 — p 2 )7rro^4(0) are the dimensionless linear con- 
ductances for the P and AP configurations, respectively. 
For the symmetric Anderson model, the Kondo effect sur- 
vives even for a finite value of polarization p (\p\ < 1), and 
g p preserves the unitary limit. As a result, the TMR is 
given by TMR = p 2 / (1—p 2 ). For the asymmetric Ander- 
son model, on the contrary, g p gets strongly suppressed 
as p increases. Then, the system exhibits a strong neg- 
ative TMR [see Fig.^c)] 0. Figure |3d) shows TMR 
as a function of e d , which shows a sharp peak around the 
symmetry point (ed ~ —U/2). The width of the peak is 
determined by how fast the Kondo effect is suppressed 
as \sd — U/2\ increases from zero, and hence depends 
strongly on the polarization p and the hybridization Tq; 
see Figs. |21 (a) and (b). Experimentally, finite tempera- 
tures would smoothen this peak. 

So far we ignored the spin-flip scattering. R ^ 
introduces a new degree of freedom in the QD, explic- 
itly lifting the level degeneracy (ed ± R). In Fig. 0(a,b) 
we plot the local DOS in the symmetric case for dif- 
ferent values of R in terms of the Kondo temperature 
k B T° K = y/rjjj2^[-Tr\ede d + U)\/2T U]. For p = 
the Kondo peak shrinks |30| and splits as expected (the 
effect can be ascribed to an external Zeeman field). For 
p = 0.25 the splitting takes place for a smaller value of 
R/Tft. As a result, Tr- in the P case is always lesser than 
T\ and TMR decreases in magnitude with increasing R 
[see Fig. Etc)]. The DOS for the asymmetric Anderson 
model are shown in Fig. 0(d,e) as a function of R/Tq as 
here we are interested in the competition between both 
splitting factors (p and R). For a nonzero value of p, 
the splitting is not symmetric unlike Fig. |3(b). This 
is again a consequence of the presence of charge fluc- 
tuations. For R — the DOS in the P alignment is 
already quenched due to the spin dependent coupling, 
which leads to a strong negative TMR. However, this ef- 
fect is washed out with increasing R and TMR tends to 
vanish [see Fig. [5(f)]. 

Conclusion. — Using a NRG method, we have shown 
that the Kondo effect in a quantum dot is not neces- 
sarily suppressed by the spin polarizations of the leads: 
for the symmetric Anderson model, where charge fluc- 
tuations are completely suppressed, the Kondo effect is 
robust even for finite polarizations. For the asymmetric 
Anderson model, the Kondo peak does split into two. We 
also reported on the TMR in the Kondo, mixed-valence, 
and empty-level regimes, being strongly affected in the 
presence of spin flip processes. The physics addressed in 



this paper is realistic and can be, in principle, tested with 
present techniques. 
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Note Added. — In the final stages of this work we be- 
came aware of a closely related work by J. Martinek et 
at, cond- mat/0304385 The difference is that they study 
only the asymmetric Anderson model, and focus rather 
on the restoring effect of an external magnetic field. 
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